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The cosmic microwave background radiation defines a preferred cosmic rest frame, and inflationary 
cosmological theories predict that the microwave background temperature fluctuations should be 
statistically isotropic in this rest frame. For observers moving with respect to the rest frame, the 

^^ , temperature fluctuations will no longer be isotropic, due to the preferred direction of motion. The 

most prominent effect is a dipole temperature variation, which has long been observed with an 
amplitude of a part in a thousand of the mean temperature. An observer's velocity with respect to 
the rest frame will also induce changes in the angular correlation function and creation of non-zero 
off-diagonal correlations between multipole moments. We calculate both of these effects, which are 

^^D ' part-in-a-thousand corrections to the rest frame power spectrum and correlation function. Both 

Cn \ should be detectable in future full-sky microwave maps from the Planck satellite. These signals will 

constrain cosmological models in which the cosmic dipole arises partly from large-scale isocurvature 

^^ ■ perturbations, as suggested by recent observations. 

^ \ PACS numbers: 98.70.Vc, 98.80.Es, 98.80.Jk 

SpH' The most prominent feature in the microwave background radiation is the large dipole modulation, at a part in 

Q . a thousand of the mean background temperature [iHal. This is generally attributed to our local peculiar velocity 
H ' with respect to the cosmic rest frame of 370 km/s [7HlC|. However, changing frames from the cosmic rest frame to a 
^ . boosted frame also induces small, distinctive changes in both the cross-power spectrum and the correlation function 
of the microwave radiation, both of which are potentially detectable in full sky maps with the angular resolution 
of the Planck satellite. Here we present a straightforward calculation of the signals and discuss their detectability, 
and note that subtle microwave background distortions are a promising route for constraining "tilted" cosmological 
'^ models where an isocurvature perturbation on the scale of the horizon contributes to the microwave dipole and to 
jy-v ' large-scale streaming motions of galaxies and galaxy clusters. 

l/~j , Similar calculations were pioneered by Challinor and van Leeuwen [11] . but they were prinrarily concerned with 

^^ ' demonstrating that the effects were small enough to be neglected when constraining cosmological parameters with 
h — , the microwave background power spectrum; they did not consider the correlation function, and did not consider 
^^ ' detectability of the signals. Buries and Rappaport [12| considered detectability of the aberration in the microwave 
^D . radiation via the shift in angular scale of acoustic peaks it introduces; this effect is related to the correlation function 
analysis we present here. 
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LORENTZ TRANSFORMATION OF THE TEMPERATURE FIELD 



A clear discussion of how a radiation field transforms under Lorentz boosts has been given by Ref. 13j, which 
clarifies some misconceptions in earlier literature, and the discussion here follows this paper. Consider a frame S' 
which is the rest frame of the microwave background, so that T'ijn!) is the temperature distribution in this frame. 
Now take frame S to be our observation frame which is boosted from the cosmic rest frame by a velocity v, with 
resulting sky temperature T(n) (we use velocity units with c = 1). Theories of cosmology predict a sky map in 
the cosmic rest frame, while we observe a sky map in our boosted frame; here we obtain the transformation of the 
correlations between multipole moments and the transformation of the correlation function between the rest frame 
and the boosted frame. 

Consider a photon with wavevector k' = k'x\! in the rest frame and wavevector k — fcn in the observation frame. 
We have the wavenumber transformation 

fc = 7(1 + vn')fc' = £»/;:' (1) 
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and dV = D~^dV' for the transformation of a volume element. We also have the aberration equation 

V • n' + t) 
v-n= (2) 

1 + vn' ^ ' 

with V a unit vector in the velocity direction. Taking derivatives of these two equations gives the transformation of 
the differential solid angle dn' = smO'dO'dcj) as dn = D^^dn'. 

In the rest frame S' , the number of photons in a volume dV' which have a wavevector fc'n' in the interval dk' dn' 
is given by 

dN = -Xt -, , \ ,, k'^dk'dV'dn' (3) 

47r3 exp(nfc'/r'(n')) - 1 

which is from the definition of a blackbody photon distribution. This same number of photons viewed in frame S will 
have propagation vector fcn in the interval dk dn and will occupy a volume dV. Rewriting Eq. ^ in terms of the 
transformed variables gives 

dN^^ ,,, , \, ,,, k^dkdVdn. (4) 

47r3 exi>{hk/DT'{n')) - 1 ^ ' 

Note that the combination k^dk dVdn is invariant, so the only change comes in the argument of the exponential. Now 
again using the blackbody definition, we identify 

T(n) = OT'(n'). (5) 

This gives the transformation between the CMB rest-frame temperature distribution T'(n') and the observation-frame 
temperature distribution T{n). From this, we now derive the transformation of observables which are commonly 
extracted from cosmological models, namely correlations of multipole moments and the angular correlation function. 

TRANSFORMATION OF POWER SPECTRA AND THE CORRELATION FUNCTION 

The (rest-frame) microwave sky temperature is commonly expressed in terms of spherical harmonics, 

T'in')^J2<nYirnin'), (6) 

Im 

where the angular power spectrum in terms of these coefficients is C/ — {a'i*^a'ij^^} . Here the angle brackets refer to 
an ensemble average over realizations of a random temperature field on the sky with the same underlying statistical 
properties. If the rest frame universe is statistically isotropic, then each moment C/ of the angular power spectrum 
is independent of m, and the average value of coefficients with different indices vanishes: (ai*m/a[„j) = if / 7^ /' or 
m ^ m' . We want the transformation law connecting the coefficients in the two frames. 
The individual aim values transform as follows: 

aim = I dnT{n)YCm{n) ^ j dn^-^^=^T'{n')Yil,{n). (7) 

Now we choose a spherical coordinate system with the z-axis aligned with the boost direction, and change integration 
variables to the rest- frame angles with n' = [9' , (f>), 

Then expanding the rest-frame temperature distribution in spherical harmonics and doing the trivial integral over 
gives the exact expression 

00 

l'=0 



pri-) = [^^7rT^] pn^) (11) 



(no sum over m) where we have defined 

/fl(„) . 2,yT^/' j^JPio:}Pr (fS) (10) 

with the abbreviation 

47r (/ + m)!^ 

for the spherical-harmonic-normalized associated Legendre functions. 

Some care must be taken in the numerical evaluation of Eq. (|10p , since the integrand is rapidly oscillating for large 
values of / — |?7i| and underflows in the integrand can compromise convergence conditions with adaptive numerical 
integrators. Direct numerical integration reveals the orthonormality relation 

E^^"i^^™. =-5,,,, (12) 

I' 

(no sum on m) which is valid for any velocity v. For small u, Eq. (|10p has the asymptotic behavior 

/;;■ ~ i + otv'). 1.^0 (13) 



«"--"<' + VItWT5)> '•^'' <"' 



'-^'' \/(2-M 2-1' ^ ' 



with all other I', I values being 0{v'^) or smaller. 

From Eq. ([9]), we have the boosted- frame products of coefficients 



I' 

Note that the statistical ensemble averaging procedure on the left side of this expression is independent of frame. 
When I2 = h, the linear perturbative expansion in Eq. (J13p is not sufficient: perturbative evaluation of Ci in the 
boosted frame requires also the 0{P'^) term for consistency, and values of / > 1/w formally require increasingly hi ghe r 
powers of v for convergence. Using the rest-frame power spectrum C/ given by the WMAP best-fit cosmology |l4|, 
direct numerical evaluation of Eq. (ITB)) with v — 0.00123 gives that Ci « 5000C( for the observed dipole, and the 
fractional corrections (Cj — C2)/C2 « 6 x lO^'^ for the quadrupole and {C[ — Ci)/Ci ranging between 10"^ and 5 x 10~^ 
for all 3 < / < 1500. These corrections to the power spectrum are too small to be observed, given the cosmic variance. 
Correlations with 1^2 ^ ^i| > 2 are at most 0{v'^) or smaller and also undetectably small as verified by direct numerical 
calculation. 

However, for the case ^2 = ^1 + 1, the linear expansions in v for /;™ arc sufficient for a consistent evaluation, yielding 



(«r..™a.™) - (QVr - C'Ml + l)^|if^ (17) 

up to corrections of O(u^). Since roughly C[ ~ l^^C2 for I up to roughly 1000 (neglecting acoustic oscillations), 
C'lj^i — C[ sa ~2C[/l and (a;*^j^ m^im) ~ ~^C!'i for large I and small to. This signal can be detected statistically, as 
shown below. 

The change in the two-point correlation function is also interesting, since distortions in the shapes of microwave 
hot and cold spots due to the boost is an effect in angle 9 space and not in multipole I space. For two sky directions 
rii and n2, the two-point correlation function is defined as 

C(ni,n2) = (T(ni)T(n2)). (18) 

In the rest frame, which we assume to be statistically isotropic, the correlation function C"(n'j^,n2) depends only on 
the angle between the two observation directions n'^ • rij . Substituting Eq. ([5]) into Eq. ([T^ yields 

C(ni, n2) ~ (1 + V • n'l + V n'2)C'{n[ ■ n'^) (19) 

up to 0{v'^) corrections. As statistical isotropy is broken by the Lorentz boost, the correlation function now depends 
on the two directions separately. The rest-frame correlation function is modified at 0{v) and varies with the angle 
between the boost direction and the observation direction. 



DETECTABILITY AND UTILITY 

For a full-sky microwave temperature map with A^pix pixels, each with Gaussian noise Upix, and a Gaussian beam of 
width (Tf,, each aim is approximately normally distributed with a variance Ci exp(— Pcr^) + w~^ [15| and uncorrelated 
with other a;,„ values, where w~^ = 47rcrp;^/A'pix is the inverse statistical weight per unit solid angle. Solving Eq. (fT7|) 
for t), averaging over / and m with signal-to-noise weighting, and propagating errors from the aim gives a standard 
error cr„ on v from a full-sky map using multipoles 2 < I < Imax of 



^l 



2 
max 



i^ i; V2(2Z + l)(2; + 3); Iv CJI Ci 



(20) 



The Planck satellite's 143 GHz channel has approximately cr(, = 3.1 arcminutes, and iVpix = 2.9x 10 with a target noise 
level of cTpix = 6.0 /iK, giving w"^ = 1.6 x 10^** /iK^. Zmax is determined by the largest I for which systematic errors 
in beam characterization do not dominate the error model for aim- For /,„ax = 2000, Eq. (I^Ul) gives cr„ = 2.5 x 10"*^. 
If the dipole is due entirely to our peculiar velocity, v = 0.00123 and Planck can detect this signal through the 
off-diagonal cross-power signal at a signal-to-noise ratio of 5. For ?niax = 2500, the signal-to- noise ratio increases to 
6; other Planck channels will provide independent estimates and further increase the signal-noise ratio. Foreground 
emission and partial sky coverage may in practice reduce somewhat the significance of a detection, although neither 
has greatly impacted measurements of the temperature power spectrum. 

Detectability of the small corrections in Eq. ([T9|) is harder to estimate, since values of the correlation function 
for similar angles are highly correlated. In general, a map contains more information about the correlation function 
than about the power spectrum, since the correlation function is significant out to separations as large as 50°. At 
a separation 9, a map has approximately {2Tr0/ab)Npiy^ pairs; for the Planck map above and, e.g., 6 = 10°, this is 
about 3.5 X 10^ pairs. Averaging over all pairs of pixels, each with Gaussian error Cpix, and propagating through the 
statistical errors on each pixel gives the standard error on C{9) as ae — CTpix -\/2C(0)/A^pairs. For a monopole and dipole- 
subtracted map, C(0) = X];(2^ + 1)C/ exp(-/2cr2)/(47r)_= 1.1 x 10** ^iK^ for the Planck beam above, so ao = 0.015 
/iK^, compared to a signal of C{9 = 10°) > 1000 fiK^ [l6|. Testing the form of Eq. (TTOt requires comparing different 
portions of the sky for a variation in the correlation function of a part in a thousand. The correlation function can 
be estimated at many different angles, with each providing a moderate signal-to-noise measurement of the difference 
in the correlation function between different sky regions. However, this estimate includes only instrumental noise, 
and does not account for cosmic variance between regions; more precise detectability estimates require evaluation of 
both the signal covariance for different angles and cosmic variance for different regions (e.g., |17|). Correlations from 
foreground emission are also a challenge for this measurement. 

Aside from being a consistency check on a fundamental cosmological property, the distinctive microwave background 
signals from a local velocity with respect to the microwave background rest frame will constrain "tilted" cosmological 
models where the dipole arises partly due to primordial superhorizon-scale isocurvature fluctuations |18l422| . Such 
models naturally explain surprising recent observations of a substantial galaxy cluster bulk flow on Hubble volume 
scales '23] and galaxy bulk flow on somewhat smaller scales 24]. 
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